Decoherence of two qubits in a non-Markovian reservoir without 
rotating-wave approximation 
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The decoherence of two initially entangled qubits in a non-Markovian reservoir has been inves- 
tigated exactly without Born Markovian approximation and rotating-wave approximation(RWA). 
The non-perturbative quantum master equation is derived and its exact solution is obtained. The 
decoherence behaviors of two qubits, initially entangled in Bell states, has been investigated in three 
different cases of parameters. The results show that the counter-rotating wave terms have great in- 
fluence on the decoherence behavior, and there are differences between the exact solution of the 
Hamiltonian with RWA and that of the exact Hamiltonian without RWA. 
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I. INTRODUCTION 

In recent years, the phenomenon , termed as 
"entanglement sudden death" (ESD), has been found 
theoretically [l|, HJ and shown experimentally [!, H|. It is 
shown that spontaneous disentanglement may take only 
a finite-time to be completed, while local decoherence 
(the normal single-atom transverse and longitudinal de- 
cay) takes an infinite timeQ]. It's quite different from 
the case of continuous variable two-atom model discussed 
by Dodd and Halliwell[5|. Thereafter, many works have 
been devoted to the related topics @, 0, H S E3l, and 
some authors extended the results in Markovian regime 
to non-Markovian case (Til. [T3|. 

The RWA, which neglecting counter rotating terms 
corresponding to the emission and absorption of virtual 
photon without energy conservation, is widely used in 
quantum optics. Generally, the RWA is justified for small 
detunings and small ratio of the atom-field coupling di- 
vided by the atomic transition frequency [l3|, H Il5l | In 
atom-field cavity systems, this ratio is typically of the 
order 10 -7 ~ 10~ 6 . Recently, cavity systems with very 
strong couplings have been discussed [16|. The ratio may 
also become order of magnitudes larger in solid state sys- 
tems, and the full Hamiltonian, including the virtual pro- 
cesses (counter- rotating terms), must be considered fl7| . 
The neglect of counter-rotating wave terms in Hamilto- 
nian strongly simplifies the mathematical treatment of 
the problem and usually give exact solution of the ap- 
proximate Hamiltonian, while the perturbative approach 
to the systems beyond the RWA usually is complicated 
and gives approximate solution[14j. Otherwise, with 
the same parameters, the approximation in Hamiltonian 
might lead to different result from the approximation of 
exact solution with exact Hamiltonian. Fortunately, non- 
perturbative master equation could be obtained by path 
integral and it is beneficial to study, exactly, the systems 



beyond the RWA [3. 

In this paper, we will focus on the influence of counter- 
rotating wave terms on the decoherence behavoir of 
two qubits in a non-markovian reservoir from the ex- 
act solution of the system beyond the RWA and Born- 
Markovian approximation. In section [H] , the reduced 
non-perturbative non-Markovian quantum master equa- 
tion of an atom in a non-Markovian reservoir is derived 
and its the exact solution is obtained by algebraic ap- 
proach of Lie superoperator. In section HTT1 the decoher- 
ence of two initially entangled atoms coupled with two 
cavities separately has been discussed. The conclusion 
will be given in section ITVl 



II. MODEL AND EXACT SOLUTION 

A. Hamiltonian and non-perturbative master 
equation 

Now we restrict our attention to two noninteracting 
two-level atoms A and B coupled individually to two en- 
vironment reservoirs [l2j]. To this aim, we first consider 
the Hamiltonian of the subsystem of a single qubit cou- 
pled to its reservoir as 



H = H a + H r 



where 
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where ujq is the atomic transition frequency between the 
ground state 1 0) and excited state 1 1) . <r z = |1)(1| — 10)(0|, 
<x_|_ = |1)(0| and <r_ = |0)(1| are pseudo-spin oper- 
ators of atom. The index k labels the field modes 
of the reservoir with frequency io k , at and a k are the 
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modes' creation and annihilation operators, and is 
the frequency-dependent coupling constant between the 
transition |1) — |0) and the field mode k. 

The reduced non-perturbative non-Markovian quan- 
tum master equation of atom could be obtained by path 
integral 

ft 



and 
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dt Pa 
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ds(Z a 







i£ (t-s) p. -i£ (s-t)\ 

^ar c jrya 

where £ a and £, ar are Liouvillian operators defined 
as 

£ p = [H a + H r ,p] 

£aP = [H a ,p] 

£arP = [H ar ,p] 

and stands for partial trace of the reservoir. 

Then we assume that the reservoir is initially in vac- 
uum states and the spectral density of the reservoir is in 
Lorentzian form[l2|, [3| 

j h = y^fffc{^( w -^fc)+<5( w +^)} 
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At 2 



2vr (w - uj ) 2 + 7 2 (6) 

where 7 respents the width of the spectral distribution 
of the reservoir modes and is related to the correlation 
time of the noise induced by the reservoir, r r = 1/7. 
The parameter A is related to the subsystem-reservoir 
coupling strength. There are two correlation functions in 
this model[20| 
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a\ comes from the rotating- wave interaction and oti from 
the counter-rotating wave interaction. 

So, we could obtain the non-perturbative master equa- 
tion of the subsystem from Eq. (J5]) 
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K , K + and K_ are superoperators 



where J , J+, J 
defined as 
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a R , a 1 and a* are real part, image part and conjugate 
of a, respectively, fit) = 1 — exp(-jt). 

From Eq.©, we find that there is a frequency shift 
caused by the interaction between the atom and the non- 
Markovian reservoir [2lj. Eq.© reveals that the contri- 
bution of counter-rotating terms is in order of A7 2 / ojq be- 
sides the rapidly oscillating terms, which is quite different 
from the result of Born and Markovian approximation of 
atom in vacuum 12211. 



B. Exact solution of master equation 



The time evolution of density operator in Eq.© could 
be obtained with algebraic approach in Ref. [23j| because 
the superoperators herein satisfy SU(2) Lie algebraic 
communication relations, i.e. 



[J_, J+] p a 
[Jo, J±] Pa 

[K.,K + ] Pa 
[K ,K ± ] Pa 
[Ki,Jj] 



-ZJoPa 
±J±Pa 

-2K aPa 
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where i,j = 0,±. By directly integrating Eq.©, the 
formal solution is obtained asflij 



p a (t) = e -r k f e ^dt(eoJ +e + J + +e-J-) 



where T is time time order operator, So = 
-t(2wo- Xja 1 ), e+ = X (7a + /(f)) /2, e_ = 
X(-ya*+f(t))/2, v Q = X ( 7 a R - /(f)), v+ = X 7 a R , 
v_ = Xf(t), T k = X {p/a R + F(t)) /2, F(t) = t - [1 - 
exp(—jt)]/j and 



/ adt = a R + ia 1 
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The exponential functions of superoperators in Eq. pip 
could be disentangled in form [3] 
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where j+, jo, j— and k+, fco, fe— satisfy the following 
equation 



A + = /!+ - /i-X 2 

A = Mo - 2p-X + 
X- = p-exp(Xo) 
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H = e for X = j and \l = v for X = fe. 

Using the results in Appendix, the exact solution of 
the master equation Eq.® is obtained 
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C. Concurrence 

In order to investigate the entanglement dynamics of 
the bipartite system, we use Wootters concurrence (24) . 
For simplicity, we assume the two subsystems having the 
same parameters. The concurrence of the whole system 
could be obtained @ 

Q = mas {0, (*,<*}, (£=*,¥) (26) 
Cl = 2e~ 2Ffc - y/pEpu) 

C2 = 2e~ 2Tk (y/pi4p4i - VP22P33) 

corresponding to the initial states of |4>) = /3|01) + 77] 10) 
and \^>) = /3|00)+7/|ll), respectively. Where /3 is real and 
< /3 < 1,7? = |r/|e^ and/? 2 + |?/| 2 = 1. The reduced joint 
density matrix of the two atoms, in the standard product 



basis <8 = {|1) = |11), |2) = |10), |3) = |01), |4) = |00)}, 
could be written as[12| 
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-217 
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here the diagonal elements are 

pu = l 2 pni0) + lmp22(0) + mlp 33 (0) - 

P22 = ippn(0) + lmp 2 2(0) + mpp 33 (0) 

P33 = ippn(0) + pmp 2 2(0) + nlp 33 (0) - 

P44 = p 2 pn(0) +pnp 2 2i0) + npp 33 (0) - 

and the nondiagonal elements are 

pu = x 2 pi 4 (0) + xyp 23 i0) 

P23 = xrp u (0) +xqp 2 s(0) 

P32 = rxp u (0) + ryp 23 i0) 

Pa = r 2 p 1A (Q) + rqp 2 3(0) 



- m 2 ( o 44 (0) 
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III. NUMERICAL RESULTS AND DISCUSSION 

In order to study the effects of non-Markovian reservoir 
on the decoherence, we assume A = IO7 in Eq.®, which 
could be realized in a hight-Q cavity ]12\. 

First, we focus on the decoherence of two qubits with 
an initial state |$). For the RWA model in Ref.dU, FigQ] 
shows that the concurrence periodically vanishes with 
a damping of its revival amplitude. For the non-RWA 
model in this paper, the decoherence of the system were 
categorized into three cases 

(A) cj > A, wo ^> 7. Fig (21 reveals that the concur- 
rence C<j> decreases exponentially to zero with very small 
amplitude oscillation and the decoherence time is about 
1 /7, which corresponds to the correlation time of ai com- 
ing from the rotating-wave terms. As for the correlation 
function 02, it has no memory effect because it's value 
averages to zero on time if j/ujq <C 1. Compared with 
Figfflfor RWA model, there is no revival of entanglement 
for non-RWA model, which exhibits that the counter- 
rotating wave terms could not be neglected in Hamilto- 
nian because the contribution of counter-rotating terms 
in master equation is in the order of A 7 2 /w 2 which influ- 
ence the decoherence behavior in long time scale. This 
characteristic will hold on for more weaker coupling con- 
stant. 
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FIG. 1: Concurrence C<s> as a function of ft and 2 with 
A = IO7 for RWA model. 








FIG. 2: Concurrence C<s> as a function of jt and 2 with 
loo = 10A, A = IO7. 



(B) luq = A > 7. From Fig{3j we could find that 
the concurrence C$ first decreases to a finite value and 
maintain it for a period of time, then periodically van- 
ishes with a damping of its revival amplitude, like that 
for RWA model. This is the result of coaction of rotat- 
ing wave process and counter-rotating process of atom 
with a non-Markovian reservoir because the contribution 
of counter-rotating wave process to the system is bigger 
than the former case and its memory effect is nozero. 

(C) ojq = 37 < A. Fig 0] exhibits that the concurrence 
C$ decreases to zero, then periodically vanishes with a 
damping of its revival amplitude, which resulted from the 
strong interaction between atom and the non-Markovian 
reservoir throuth virtual photon process. 

From the discussion above, we find that the decoher- 
ence behavior of C$ is symmetric with (3 2 because of the 
symmetry of initial state |$). And the revival of entan- 
glement could be found only with big ratio of coupling 
constant to the atom transition frequency. For no-RWA 
model, the decoherence behaviors are richer than that for 




FIG. 3: Concurrence C$ as a function of jt and with 

LUO — X — IO7. 




FIG. 4: Concurrence C<s> as a function of "ft and f3 2 with 
A = IO7, wo = 37. 



RWA model. 

Then, we focus on the decoherence of two qubits with 
initial state of |#). For the RWA model in Ref.0, Figd] 
shows that the entanglement represented by GV has a 
similar behavior to C$ for (3 2 > 1/2 in Fig[TJ In con- 
trast, for 1 < 1/2, there is ESD because vanishes 
permanently after a finite time, similar to the Markovian 
case0- Second, revival of entanglement appears after 
periods of times when disentanglement is complete. For 
the non-RWA model, the decoherence of the system were 
also categorized into three cases 

(A) ojo > A, lj » 7. Fig [6] reveals that the concurrence 
GV decreases exponentially to zero for (3 2 > 1/2, while 
CV vanishes permanently after a finite time for /3 2 < 1/2. 

(B) loo = A > 7. FigJT] shows that the entanglement 
represented by CV has a similar behavior to C$ for f3 2 > 
1/2 in FigO In contrast, for (3 2 < 1/2, the concurrence 
first decreases to zero and maintain it for a long period 
of time, then revives with very small amplitude before 
vanishes permanently. 
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FIG. 5: Concurrence C* as a function of -yt and 1 with 
A = IO7 for RWA model. 




FIG. 6: Concurrence CV as a function of ft and 2 with 
loo = 10A, A = IO7. 




FIG. 7: Concurrence CV as a function of yt and 2 with 
lo = A = IO7. 



(C) wo = 37 < A. The evolution dynamics of concur- 
rence C^i is almost the same as that in Fig0J Unlike 
the two cases above, the evolution behavior of concur- 
rence Cijf becomes symmetric, like that of C$, because 
of the strong interaction of atom with reservoir through 
the emission and absorption of virtual photon. 

The above characteristics of also verify that the 
revival of entanglement could be found only with big ratio 
of coupling constant to the atom transition frequency. 

IV. CONCLUSION 

The reduced non-perturbative non-Markovian quan- 
tum master equation of atom in non-Markovian reservoir 
has been derived and its exact solution is obtained by al- 
gebraic approach of Lie superoperator. The decoherence 
of two initially entangled atoms, coupled with two vac- 
uum cavities separately, has been discussed. 

The results show that the decoherence behavior of two 
qubits in a non-markovian reservoir is dependent on the 
ratio of the coupling strength to atomic frquency. First, 
with the increasing of coupling strength, the decoher- 
ence behavior becomes more and more complicated, and 
there is revival of entanglement after a period of time of 
disentanglement, due to the strong interaction and the 
memory effect of the non-markovian reservoir. Second, 
the strong coupling and the counter-rotating wave inter- 
action could make the decoherence of even parity Bell 
entanglement state become symmetrical. Third, with 
no-RWA model, one could find much more decoherence 
behaviors than that of RWA model under different con- 
ditions of system parameters. The case for finite temper- 
ature could also be obtained by the method in this paper 
and will be published elsewhere. 
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APPENDIX A 

Defining two superoperators as 

J = e 3+J+ e hJo e j-J- ( A1 ) 

K = e k+K +e k K 0e k-K- ( A2 ) 

and using the following relations 

e j + J +p = p + j+a +P <7 + (A3) 

e j - J -p = p + j-a-pa- (A4) 

e k+K +p = p + k+a+pa- (A5) 

e k - K -p = p + k-a-po + (A6) 
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(A7) 
(A8) 



one could get 



pu pio \ ( Pn {e j °/ 2 + e la,2 j+i-)pw + e Jo/2 j+p i 



Poi Poo I V e Jo/2 p i + e 3a,2 j-Pw Poo 



(A9) 
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Kl PoiAx>J~(v poi e-W2 poo + e -W2 fc _ pil I 
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